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ABSTRACT 

We present Chandra X-ray observations of the nearby Seyfert 1.5 galaxy 
NGC 4151. The images show the extended soft X-ray emission on the several 
hundreds of pc scale with better sensitivity than previously obtained. We show 
that the hard X-ray component (> 2 keV) is spatially unresolved. The spec- 
trum of the unresolved nuclear source may be described by a heavily absorbed 
{Nh — 3 X 10^^ cm~^ ), hard power-law ( F ~ 0.3) plus soft emission from either 
a power-law (F ~ 2.6) or a thermal {kT ~ 0.6 keV) component. The flux of the 
high energy component has decreased from that observed by ASCA in 1993 and 
the spectrum is much harder (F ~ 0.3 in 2000 versus F ~ 1.5 in 1993). The large 
difference between the soft and hard spectral shapes does not favor the partial 
covering or scattering model of the "soft excess" . Instead, it is hkely that the 
hard and soft nuclear components represent intrinsically different X-ray sources. 
The stronger nuclear emission lines of those seen by the Chandra HETGS spec- 
trum are detected. Spectra of the extended emission to almost 1 kpc NE and 
SW of the nucleus have also been obtained. The spectra of these regions may be 
described by either thermal bremsstrahlung [kT ~ 0.4 — 0.7 keV) or power-law 
(F ~ 2.5 — 3.2) continua plus 3 emission hues. There is an excellent correlation 
between the extended X-ray and [O lll]A5007 line emissions. We discuss the 
nature of the extended X-ray emission. Because there is no extended electron- 
scattered hard X-ray emission, an upper limit to the electron scattering column 
can be obtained. This upper limit is much too low for the soft X-rays to be elec- 
tron scattered nuclear radiation, unless the nucleus radiates soft X-rays much 
more strongly towards the extended regions than towards Earth, a situation we 
consider unlikely. We favor a picture in which the extended X-ray emitting gas is 
heated in situ by the nuclear radiation. Some of the X-rays may originate from 
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a hot phase which confines the warm ionized gas seen optically, although X-ray 
emission produced via photoionization by the nucleus is also likely. A faint, prob- 
ably background, compact X-ray source lies ~ 2f2 from the nucleus to the SW, 
approximately along an extension of the extended SW X-ray and [OIII] emission. 

Subject headings: galaxies: active - galaxies: individual (NGC 4151) - galaxies: 
nuclei - galaxies: ISM - galaxies: Seyfert - X-rays: galaxies 

1. Introduction 

Because of its proximity (13.2 Mpc for Hq = 75 km s~^ Mpc^^, 1" = 64 pc), the Seyfert 
1.5 galaxy (Osterbrock & Koski 1976) NGC 4151 is one of the best studied among the class. 
X-ray observations show that NGC 4151 is a modest luminosity (-^2-10^6^ 2 — 20 x 10^^ 
gj,gg-i ^ active galactic nucleus (AGN) with a highly variable and hard continuum. There is 
evidence for a spectral-fiux correlation above 2 keV (Perola et al. 1986; Yaqoob & Warwick 
1991; Yaqoob et al. 1993). The soft X-ray continuum of NGC 4151 is characterized by a 
prominent excess of X-ray emission above an extrapolation of the absorbed, higher energy 
power-law component to lower energies. This "soft excess" has been described in terms 
of partial covering, a warm absorber, scattered nuclear radiation and additional spectral 
components (e.g. Holt et al. 1980; Yaqoob & Warwick 1991; Weaver et al. 1994a, b). 

Einstein and i? 05^1 T observations identified extended soft X-ray emission that appears 
to be associated with ionized gas visible optically (Elvis, Briel & Henry 1983; Morse et 
al. 1995). This extended emission must account for at least part of the "soft excess". Ogle 
et al. (2000) have recently obtained a 48 ks observation of NGC 4151 with the High-Energy 
Transmission Grating Spectrometer (HETGS) on board the Chandra X-ray Observatory. 
Their observation provided the first high resolution X-ray spectroscopy of NGC 4151 and 
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also a zero-order image. Strong, narrow X-ray emission lines are seen in the spectrum. 
Evidence of both photoionized and coUisionally ionized gas was found. However, their 
zero-order image has only modest signal-to-noise due to the limited transmission of the 
gratings. 

We thought it worthwhile to obtain an exposure without the gratings in place in order 
to obtain a sensitive image of the extended emission. Given the dependence of the HETGS 
zero-order transmission on energy, our ~ 30 ks direct image has higher sensitivity than the 
zero-order image by Ogle et al. (2000) for all energies below ~ 6 keV. 

In this paper, we present a Chandra Advanced CCD Imaging Spectrometer (ACIS) 
observation of NGC 4151. Observations are described in Section 2. Data analysis and 
results are presented in Section 3, where we present the best X-ray image of NGC 4151 
to date and low resolution spectra of both the nucleus and extended regions. The X-ray 
image is also compared with optical and radio images. Implications of the observation are 
discussed in Section 4. The detection of an off nuclear source about 2^2 SW of the nucleus 
of NGC 4151 is described in the Appendix. 

2. Observations and Effects of Pile Up 

The observations of NGC 4151 were obtained in two phases. First, to evaluate the 
effects of pile up^ in the bright inner regions of NGC 4151, we took a preliminary 4.3 ks 
observation on Dec. 4 1999 in so called "alternating mode", in which single exposures 
with a 0.1 s frame time were alternated with two exposures with a 0.4 s frame time. The 

^For bright sources, two or more photons may fall in the same pixel during one frame time 
interval and the photons are then detected as a single event. This "pile up" effect results in 
an underestimated count rate and an incorrect spectrum. 
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observation showed that the nucleus is significantly piled up in the 0.4 s, but not in the 0.1 s, 
frame time observations. We also inferred that the bright inner regions of the galaxy would 
be strongly affected by pile up in the planned long exposure, which would use the standard 
3.2 s frame time. We thus decided to obtain another, longer observation in "alternating 
mode", in addition to the long observation, to mitigate the effect of pile up in the inner 
regions. Therefore, the second phase of observations consisted of a ~ 26 ks observation 
with 3.2 s frame time and a ~ 7 ks observation in "alternating mode", including both 0.1 s 
and 0.4 s exposures. Observations are summarized in Table 1. In the ~ 26 ks observation, 
CCDs 12, 13, SI, S2, S3 and S4 were read out, though all of the detected X-ray emission 
from NGC 4151 is on S3. In both "alternating mode " exposures, only chip S3 was read 
out. 

The fraction of bad grades (i.e. the ratio of the number of counts with ASCA grades 
1+5+7 [defined in the the Chandra Proposers' Observatory Guide, section 6.3] to the total 
number of counts) was found to be a good indicator of the degree of pile up in any given 
pixel through comparison with the event rates in the different frame time observations. 
Regions with a large fraction (nominal value > 10%) of bad grades were considered to be 
piled-up. As already found in the preliminary observation, the 0.1 s frame time data were 
not piled-up, with less than 6% of the grades in the brightest pixel being bad. The 0.4 s 
and 3.2 s frame time data were found to suffer from pile-up out to 1'.'3 and 3" from the 
nucleus, respectively. 

All of the analysis was initially done with the "old" response functions ( "fefs" ) , and 
data were reduced and analysed with CIAO vl.1.5 and XSPEC vll.0.1. The observations 
were screened for high background count rates and aspect errors by following the Chandra 
Science Threads. After we received the referee's report, new versions of response function 
(CALDB 2.7) were available. We reanalysed some of the data with these new response 
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functions and with the latest version of CIAO - v2.1.2. In general, the change in derived 
parameters (such as absorbing column, photon index, bremsstrahlung temperature and 
emission line properties) were small. We did not reanalyse the 0.1 s frame time nuclear 
observations with the new response functions because our results agree well with the 
concurrent HETGS observations of the nucleus by Ogle et al. (2000), which have much 
higher spectral resolution and are thus superior for detection of lines. We have, however, 
added 0.2 keV to the line energies below 1 keV from the 0.1 s frame time observation since 
there is a systematic error of this magnitude when using the old response functions. The 
results presented for the 0.4s frame time nuclear observations and the 3.2s frame time 
observations of the extended emission employed the new response functions. 

3. Data Analysis and Results 

3.1. Morphology of X-ray emission 

The soft X-ray (0.3-2.5 keV) images derived from the 0.1 s, 0.4 s and 3.2 s frame 
time exposures are shown in Fig. 1. The "hole" at the position of the nucleus in the 3.2 s 
frame time image (Fig. Ic) is due to pile-up. The X-ray emission of NGC 4151 comprises 
a bright unresolved nucleus and resolved extended regions. The south-west part of the 
extended emission extends as far as 14" ( ~ 900 pc) from the nucleus along PA ~ 233°. 
The north-east part of the extended emission extends to 11" (~ 700 pc) from the nucleus 
along PA ^ 75°. The structure of the extended emission appears to be knotty rather than 
smooth. A moderately bright, compact knot is clearly visible 6" (~ 380 pc) south-west of 
the nucleus (Fig. Ic). A contour map representation of the extended emission is shown in 
Fig. 2. Our images are in good agreement with the HETGS zero-order image obtained by 
Ogle et al. (2000), but are more sensitive. 
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The 2-9 keV emission is not resolved in our observations, as can be seen by comparing 
the radial profile of the observations in this band with the point spread function (PSF). 
In Fig. 3a, the 5.5 keV model PSF without aspect error is plotted with the observed 
profiles of the 2-9 keV emission. The 0.4-2 keV emission is, on the other hand, resolved 
along the northeast-southwest direction, but is not (or is marginally) resolved along the 
northwest-southeast direction (Figs. 3b, c, d). 

3.2. X-ray Spectra of the Nuclear Region 

We extracted spectra of the nucleus from a circular region with radius 3" in the 0.1 s 
and 0.4 s frame time observations. The total number of counts in the spectra are 2033 
and 10690 in the 0.1 s and 0.4 s frame times, respectively. For both observations, the 
background was taken from a region between two ellipses centered on the nucleus. The 
large ellipse had a semi-major axis of 17'.'5 along PA 27° (the column direction of CCD read 
out) and the smaller ellipse had a semi-major axis of 12'.'5 in the direction of PA = 52°. This 
rather oddly-shaped background region allowed all the background events to be extracted 
from the same CCD node. 

The extracted spectra were grouped to a minimum of 20 and 30 counts per bin in the 
0.1 s and 0.4 s frame time data, respectively, to allow use of statistics. The spectral 
resolution of ACIS-S is ~ 0.1 keV at 1 keV. For the bright nuclear region, the spectral 
resolution was not significantly degraded by this grouping. Because of the poor calibration 
at low energies and high background at high energies, channels below 0.3 keV and above 9 
keV were ignored when modeling the spectra. 

The emission from the brightest part of NGC 4151 originates in the nucleus and the 
unresolved emission line regions. At least two components are needed to describe the 
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continuum of the spectra from this region: a hard power law absorbed by an equivalent 
hydrogen column density of Nh = (2 — 3) x 10^^ cm~^ and a soft component that can be 
represented by either an absorbed power law or a thermal bremsstrahlung model. The 
0.1 s frame time spectra were fitted with these two component continuum models (Table 2) 
plus emission lines (Table 3). The model 2-9 keV flux is 4.8 x 10~^^ ergcm~^s~^ and the 
photon index of the hard component F ~ 0.32. The soft X-ray spectrum below 2 keV can 
be described as a power law with photon index F ~ 2.6 (Table 2, Fig. 4a). An alternative 
description in terms of thermal bremsstrahlung gives a temperature ~ 0.57 keV (Fig 4b). 
However, the emission cannot be described by solar abundance, coUisionally ionized hot 
gas. The best fit double power law parameters agree with those obtained in a concurrent 
HETGS observation (Ogle et al. 2000). Previous ASCA observations (Weaver et al. 1994b) 
gave a 2-10 keV flux of 11 — 22 x 10"^^ ergcm~^s~^ , a photon index of F ~ 1.5 and an 
absorbing column of 2 — 5 x 10^^ cm~^ . Weaver et al.'s absorbing column is about the 
same as found in our observation and that of Ogle et al. (2000), while the photon index and 
flux from ASCA are signiflcantly larger, indicating variability of the hard component. Our 
observations suggest that the hard component of NGC 4151 is at its lowest state since 1984 
(EX05'^r observations yielded a low 2-10 keV flux of 3.6 x 10~^^ ergcm~^s~^ , [Pounds et 
al. 1986]). The very hard power law index found in the present observations with Chandra 
is unusual for a Seyfert 1 galaxy (Mushotzky 1984). 

Lines were identifled by adding narrow Gaussian features to the continuum model 
with the goal of reducing the value of the flt. For all cases except the combined O Vll 
triplet, the intrinsic line width was set to zero. Both 0.1 s and 0.4 s frame time data were 
used. Though piled-up, the continuum of the 0.4 s frame time data between 0.4 to 7 keV 
can still be described by a double power law model reasonably well (reduced — 180.8 
for 179 degrees of freedom), with a rather hard photon index (F = — 0.14lo;o7) above 2 
keV (Fig. 4c). This value of F is, however, too low as a result of the pile up. The low 
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spectral resolution of the ACIS instrument does not allow us to resolve all the lines seen 
in the grating observation (Ogle et al. 2000), but the stronger hues are detected. The best 
fits to the spectra are shown in Fig. 4, while the lines identified are listed in Tables 3 and 
4. The errors represent the 90% confidence range for a single interesting parameter. One 
feature in the 0.4 s frame time spectrum, which probably results from the rapid change of 
effective area with energy near 2 keV, is also listed in Table 4 and labeled as "artefact" . 
The line fluxes from the 0.1 s frame time image agree well with those found in the HETGS 
observation (Table 3). The hne fluxes found from the 0.4 s frame time data tend to be 
lower than those obtained from the grating observation (see Table 4) , probably as a result 
of pile up. 

3.3. X-ray Spectra of the Extended Emission 

For the extended emission, we used the 3.2 s frame-time observation and extracted 
spectra from non-piled-up regions more than 3" from the nucleus to both NE and SW. The 
spectrum of the NE region was extracted from a 12^'7 x 8''6 rectangle centered 7'! 3 northeast 
of the nucleus {a = 12^* 10™33*.0, 6 = 39° 24' 26'.'7) with the longer side in position angle 
144°. The spectrum of the SW region was extracted from a 12'/3x7'.'7 rectangle centered 
9'/2 southwest of the nucleus {a = 12'* 10"^ 31^8, 6 = 39°24'ir.'9 ) with the longer side 
in position angle 234° (Fig. 5). The total number of counts in the NE and SW regions 
are 1120 and 1776, respectively. The background region was the same as in section 3.2. 
The data were grouped to a minimum of 20 counts per bin in each region. The spectral 
resolution above 1 keV was significantly reduced by this grouping, especially for the NE 
region. As before, channels below 0.3 and above 9 keV were ignored in spectral modelling. 

The extended soft X-ray continuum emission can be modeled by a single bremsstrahlung 
or power law model with photo-electric absorption (Table 5). The emission observed above 
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2 keV is the unresolved nuclear emission which has been spread out by the mirror PSF 
(see section 3.1 and Fig. 3a), and are hence modeled with an absorbed power law with 
photon index and absorbing column density identical to the nuclear value (section 3.2). 
The absorbing column densities for the extended regions are not well constrained with the 
present data and were "frozen" to the Galactic value (Nn.cai = 2.19 x 10^° cm~^ , Murphy 
et al. 1996). 

It would be interesting to examine the spectral index as a function of distance from the 
nucleus. However, the low count rate in the extended region makes it hard to obtain such 
spectra. To obtain a more qualitative picture, the "hardness ratio" of the soft X-ray band 
(defined as the ratio of background subtracted count rates in the 0.8-2.0 keV and 0.4-0.8 
keV bands) as a function of distance from the nucleus has been obtained for the SW region 
from the 3.2 s frame time observation (Fig. 6). No significant change of hardness ratio is 
seen between 4" and 14" from the nucleus. 

Spectral lines were identified with the same technique as for the nucleus, and the 
results are shown in Table 6. He-hke O and Ne triplets are seen in both the SW and NE 
regions. The 0.72-0.74 keV feature observed in both regions could be a blend of lines of 
H-like O and radiative recombination continuum features of O VII. The observed spectra 
are compared with the continuum plus emission-line models in Fig. 7. 

3.4. Compcirison with observations at other wavelengths 

An [O lll]A5007 contour map ( Perez- Fournon & Wilson 1990) is overlaid on the 3.2 s 
frame time 0.3-9 keV image in Fig 8. Extended X-ray emission is well correlated with the 
forbidden line emission to both NE and SW of the nucleus. The X-ray knot 6" (^ 380 pc) 
to the SW of nucleus is also seen in the [O ill] image, but with a ~ 1" displacement. 



- 11 - 



In Fig. 9, a radio contour map from MERLIN observations (Mundell et al. 1995) is 
overlaid on our 0.1 s frame time Chandra image. The radio jet does not align with the soft 
X-ray extension. The angle between the two directions is ~ 12°. 

4. Discussion 

4.1. Pctrtial covering and Scattering 

One of the early suggestions for the origin of "soft excess" is the partial covering model 
(Holt et al. 1980), in which the continuum from the nucleus is covered by a "clumpy" 
absorber. Previous studies show that a covering factor of 0.7 to 0.97 can account for all 
or a significant part of the "soft excess" (e.g. Weaver et al. 1994b, Warwick et al. 1995). 
However, in the two concurrent observations presented by Ogle et al. (2000) and this paper, 
the photon indices for the hard and soft component were found to be ~ 0.3 and ~ 2.6, 
respectively. The 0.1-2 keV flux from our observation is ~ 4 x 10~^^ ergcm~^s~^ , which is 
close to the lowest value seen in the ASCA observations of ~ 5 x 10~^^ ergcm"^ s~^ (Weaver 
et al. 1994b). This shows no significant variabihty in the soft component, in contrast to 
the large variability in the hard component. These results makes a partial covering model 
unlikely. The same argument applies to the notion that the soft emission is Thomson 
scattered nuclear radiation, as already noted by Ogle et al. (2000). Instead, it is most likely 
that the soft and hard components are intrinsically different X-ray sources. 

The Thomson scattering model for the extended soft emission can be evaluated using 
the soft- and hard-band profiles (Fig. 3). The hard-band profile is consistent with the mirror 
PSF, indicating there is no detectable intrinsically extended hard emission. This allows an 
upper limit to the fractional brightness of any hard nuclear emission which is Thomson 
scattered by the extended gas to be obtained. Because Thomson scattering is wavelength 
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independent, this upper limit can be compared with the brightness of the extended, soft 
emission, expressed as a fraction of the soft nuclear emission. Using the 0.1 s frame time 
observations in Fig. 3a, the 2-9 keV brightness within 1" of the nucleus is 0.12 cts pixel"^ 
s"^ and the intrinsic 2-9 keV brightness between 4" and 5" from the nucleus is < 3.9 x 10~^ 
cts pixel"^ s""^, which represent an upper limit to the Thomson scattered hard radiation. 
The corresponding numbers for the 0.4-2 keV emission are 0.13 cts pixel"^ s~^ for the 
nucleus and 2.6 x 10~^ cts pixel"^ s~^ between 4" and 5" from the nucleus in the SW region 
(Fig. 3c). Defining R =(brightness between 4" to 5" from the nucleus)/ (brightness within 
1" of the nucleus), we have -R2-9 keV < 3.2 x 10~^ and Rqa-2 keV = 2 x 10~^. The argument 
indicates that the extended soft emission is not Thomson scattered nuclear radiation as long 
as the intensity of the nuclear radiation radiated towards the putative scattering region is 
equal to that radiated towards Earth. The nuclear 2-9 keV radiation radiated towards the 
extended region may be stronger than that radiated in our direction because the latter is 
absorbed. Such an effect decreases R2-9 keV and increases the discrepancy between R2-9 keV 
and RoA-2 keV- If the soft nuclear emission is absorbed by more than the low column we 
obtained by fitting the spectrum, it is possible that more soft radiation is radiated from 
the nucleus towards the extended emission than towards us. This effect would reduce 
RoA-2 keV a-nd could potentially remove the discrepancy between Rqa-2 keV a-nd R2-9 keV- 
While such cannot be ruled out, we consider it unlikely that the extended X-ray emission is 
Thomson-scattered nuclear radiation. 

4.2. A Two Phase Model for the Extended X-ray Emission 

The close association between the extended soft X-ray emission and the optical 
hne emission in NGC 4151 suggests that the X-rays may arise from hot gas in pressure 
equilibrium with the extended narrow line region (ENLR) (Heckman & Balick 1983, Morse 
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et al. 1995). Gas of 10^ — 10^ K is needed to confine the cooler gas clouds (Krolik, McKee 
& Tarter 1981, hereafter KMT; Krolik k Vrtilek 1984). Grating observations with the 
Chandra HETGS show clear signatures of both photoionization and coUisional ionization in 
the X-ray emitting gas in the narrow line region. This indicates at least two phases exist 
in the plasma, one being hot (T ~ 10^ K) and coUisionally ionized and the other warm 
(T > 10^ K) and photoionized (Ogle et al. 2000). 

Penston et al. (1990) derived emissivity- weighted estimates for the electron density 
and electron temperature in the SW optical ENLR of — 220 cm~^ and Tg ~ 14130 
K, corresponding to a pressure of UeTg — 3.1 x 10^ cm~^K. If the warm and hot clouds 
are in a rough pressure equilibrium, the density of the hot phase gas at a temperature of 
Th — 4.8 X 10^ K (the bremsstrahlung temperature in the SW regions) should be Uh ^ 0.6 
cm~^ . Assuming the hot gas in the observed SW region uniformly fills a cylinder with a 
radius of 246 pc (~ 3'.'85) and a height of 787 pc (~ 12'/3 ), the number density inferred 
from the bremsstrahlung description of the continuum of the SW radiation (see Table 5) 
is found to be Uh — 0.5 cm~^ . This is close to that needed for pressure balance. Using 
the two power- law spectrum obtained from the 0.1 s frame time observation (section 
3.2) and adopting a spectral turn-over near 100 keV (Maisack et al. 1993), we obtain 
a photoionizing luminosity Li^s ev-ioo keV = 6.2 x 10^^ ergcm~^s~^ . The ionization 
parameter (defined as 5 = L/iirckTrir'^, KMT) for the hot gas measured at the center of 
the SW region (r = 589 pc) is then 0.19. This combination of temperature and ionization 
parameter lies in the general area of the equilibrium curves of T versus S calculated by 
KMT, although the equilibrium would be unstable. We note that the temperature derived 
from the bremsstrahlung model is much lower than the Compton temperature ( T > 10^ K, 
see KMT). Given our low spectral resolution, we cannot tell whether the emission lines we 
see from the extended region are coUisionally ionized or photoionized. Ogle et al. (2000) 
identified photoionized N vil, O Vll, O Vlll, Ne ix and Ne x in the nuclear region from 
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the narrowness of their radiative recombination continua. We see lines (some tentatively) 
of O VII, Ovill, Ne IX and Ne X from the extended regions (Table 6). If these species are 
created by photoionization, the ionizing parameter lies in the range 1.6 < log ^ < 2.2, where 
^ = L/nr"^ (Kallman & McCray 1982). Using the ionization luminosity given above, and a 
typical radius of 500 pc, these values of ^ imply 0.2 < n < 0.7 cm~^ . The temperature 
of this gas would be ~ 10^'^ K according to Kallman & McCray (1982). Higher spectral 
resolution observations are needed to determine whether these lines are photoionized in 
such a warm gas or coUisionally ionized in a hotter one. 

5. Conclusion 

We have obtained high sensitivity X-ray imaging-spectroscopy of NGC 4151 with the 
Chandra X-ray Observatory . The soft X-ray emission (below 2 keV) is spatially resolved 
and extends as far as ~ 900 pc to the SW and ~ 700 pc to the NE of the nucleus. There 
is a close correlation between the extended soft X-ray and [O lll]A5007 emissions on the 
hundreds of pc scale. The 0.4-9 keV spectra of the nucleus and the extended regions to 
the SW and NE of the nucleus have been obtained. The X-ray emission above 2 keV is 
spatially unresolved and is well described by an absorbed power law. The spectrum above 2 
keV is harder and the flux lower than in observations made with ASCA in 1993 (Weaver et 
al. 1994b). The spectrum below 2 keV, on the other hand, is softer than found previously. 
The large difference between the soft and hard spectral shapes does not favor the partial 
covering or scattering model of the "soft excess" . Instead, it is likely that the hard and soft 
nuclear components represent intrinsically different X-ray sources. We have shown that a 
model in which the extended emission is electron-scattered nuclear radiation is extremely 
implausible. With the temperature (~ 4.8 x 10^ K) inferred from a bremsstrahlung 
interpretation of the X-ray spectrum, this hot gas is close to pressure equilibrium with 
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the gas which radiates the optical hne emission. While the extended emission is certainly 
thermal, our observations do not permit a clear distinction between a photoionized and 
coUisionally ionized gas. It is likely that both components are present, given the findings of 
Ogle et al. (2000) for the nuclear region. 

We thank Andy J. Young, Patrick L. Shopbell and David A. Smith for help with 
using CIAO, XSPEC and WIP in analyzing the data. We are also grateful to Carole G. 
Mundell for providing the MERLIN map in the numerical form. This research was partially 

supported by NASA grant NAG 81027, by HST grant GO 7275, and by the Graduate 
School of the University of Maryland. 

APPENDIX: An Off Nuclear X-ray Source 

An unresolved X-ray source (q;(2000.0) = 12'* 10"*22^37, 5(2000.0) = 39° 23' ir/0) ~2;2 
SW of the nucleus (P.A. ~ 241°) of NGC 4151 has been detected in the 3.2 s frame time 
observation. The total counts within a 3" radius centered on the point source is 249. We 
obtained a spectrum of the source (Fig. 10). The spectrum can be well described by a 
power law absorbed by the Galactic column density (Nn.oai ^ 2.19 x 10^° cm~^ , Murphy et 
al. 1996), with = 27.2 for 25 degrees of freedom. The best fit photon index F = l.ltol 
and the model fiux for the 0.3-5 keV energy range is 6.8 x 10~^^ ergcm"^ s~^ , corresponding 
to a luminosity of 1.4 x 10^^ ergs~^ , if the source is associated with NGC 4151. This 
source lies approximately along the extension of the SW extended emission to larger radii. 
There is a weak optical object nearby at the limit of the DSS with an estimated POSS-E 
magnitude of ~ 20. It is unclear whether this X-ray source is physically associated with 
NGC 4151. Similar coincidences have been seen in Mkn 3, Pictor A and NGC 4258, each of 
which has a faint, compact X-ray source close to the extension of the radio jets (Morse et 
al. 1995; Wilson, Young & Shopbell 2001; Wilson, Yang & Cecil 2001). 
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Table 1. Observation Summary 



Obs ID 


Sequence # 


Observation date 


Frame time 


Exposure^ 


347 


700019 


12/04/99 11:29:54 


0.1 sec 


296.8 sec 


347 


700019 


12/04/99 11:29:54 


0.4 sec 


2374.3 sec 


348 


700020 


03/15/00 10:51:40 


3.2 sec 


26.2 ksec 


372 


700198 


03/06/00 13:27:39 


0.1 sec 


771.1 sec 


372 


700198 


03/06/00 13:27:39 


0.4 sec 


6168.8 sec 



''Total good time with dead time correction. 



Table 2. Spectral Models of the Nucleus Prom the 0.1 s Frame Time Observations 



ModeP 


Nh (Soft Comp.) 
[xl020cm-2] 






iVij (Hard Comp.) 
[xl022cm-2] 


r2 




/ d.o.f. 


PL+PL+LINES 

BR+PL+LINES 


2.01°;^ 




(7.8tO;^) X 10-4 
{2.9+J,l) X 10-3 


o 1+0.4 
"'•■'--0.3 


n 09+0.05 
'-'•"^^-0.12 

U.OZ_g 21 


(i.4iE5:J) X 10-3 
(1.410:3) X 10-3 


66.3/65 
70.4/69 



"^Model abbreviations are: PL = power law, BR = bremsstrahlung, LINES = emission lines (see Table 3 for line parameters). 

''Unit: kcV 

'^Modcl normalization. For the bremsstrahlung model XBrcm = 3.02 x 10^^^/ netiidV j {'i'^D'^), where n,, is the electron density, n/ 
is the ion density and D is the distance to the source (all in cgs units). For the power law model Kpi^ = phcm" -^s-^keV-^ at 1 keV. 



Table 3. Emission Lines Prom the Nucleus: 0.1 s Frame Time Observations 



ModeP 


Energy 


Line 


Observed Energy 






Flux<i 




[keV] 




[keV] 




[eV] 





I 


6.40 


Fe I Ka 


f. 07+0.04 

""J' -0.06 




X 10" 


-4 


182 


1.8 X 10" 


-4 


I 


0.87, 0.90, 0.91 & 0.92 


VIII RRC, Ne ix triplet 


n 01+0 01 
u.yi_o.o2 


(8.7tl1) 


X 10" 


-5 


78.0 


7.4 X 10" 


-5 


I 


0.74 & 0.77 & 0.82 


VII RRC & Vlli Ly/3 Ly7 


n 7(5+0.06 
^- '"-0.05 




X 10" 


-5 


21.7 


8.1 X 10" 


-5 


P 


0.561, 0.568 & 0.574 


VII triplet 


^•3°-0.04 




X 10" 


-4 


71.0 


4.3 X 10" 


-4 


II 


6.40 


Fe 1 Ka 


6 36+0 04 




X 10" 


-4 


181 


1.8 X 10- 


-4 


II 


0.87, 0.90, 0.91 & 0.92 


VIII RRC, Nc IX triplet 


qi+o oi 

^•^-'^-0.01 




X 10- 


-4 


105. 


7.4 X 10- 


-5 


II 


0.74 & 0.77 & 0.82 


VII RRC & VIII Ly/? Ly-f 


76+0 01 

^- '"-0.03 




X 10" 


-5 


50.7 


8.1 X 10- 


-5 


II 


0.65 


VIII Lya 


62+" °2 


(l.6tS]-J) 


X 10- 


-4 


55.2 


1.0 X 10- 


-4 


II 


0.561, 0.568 & 0.574 


VII triplet 


t^s+o-o^ 


(l.7irr) 


X 10" 


-4 


37.4 


4.3 X 10" 


-4 



'^Model I: bremsstrahlung + power law + emission lines; Model II: two power laws + emission lines. 

= total phcm~^ in the line. 
'^Equivalent Width 

•^Total phcm-2 s'^ in the line from HETGS observations by Ogle et al. (2000) 

•^Modelled as a single line with width 0.05±0.03 keV; this width is probably not real in view of present calibration uncertainties. 



Table 4. Emission Lines Prom the Nucleus: 0.4 s Frame Time Observations 



Energy 


Line 


Observed energy 




a 








Flux'^ 




[keV] 




[keV] 










[eV] 






6.40 


Fe I Ka 


fi 4-3+0.09 


(5.9t|°o) 


X 


10- 


-5 


92.6 


1.8 X 10- 


-4 


2.25 


Artefact 












- 


- 




1.73, 1.74, 1.84 & 1.88 


Si I Ka, Mg XII Ly/3, Si xiii f, r 


■'-•O-'--0.02 


_ _i_n fi\ 


X 


10" 


c 
-0 


146. 


4.0 X 10- 


c 
-0 


1.33, 1.35 & 1.36 


Mg XI f, r, Ne x RRC 


1 qo+0.03 

1.00_Q Q3 


(7.lttS) 


X 


10- 


-6 


25.7 


2.9 X 10- 


-5 


1.02 


Ne x Lya 


^•^'-0.03 


(2.0;;,;^) 


X 


10- 


-5 


39.9 


2.0 X 10- 


-5 


0.90, 0.91 & 0.92 


Nc IX triplet 


0.921°:- 


(6.9;}j) 


X 


10- 


-5 


100. 


6.2 X 10- 


-5 


0.74, 0.77 & 0.82 


VII RRC, VIII Ly/3, Ly7 


0.8lt°f3 




X 


10- 


-5 


54.5 


8.1 X 10- 


-5 


0.65 & 0.66 


VIII Lya, N vii RRC 


'^•"'-0.02 




X 


10- 


-5 


60.8 


1.2 X 10- 


-4 


0.561, 0.568 & 0.574 


Vll triplet 


56+° °2 


(3.51^:^) 


X 


10- 


-4 


166 


4.3 X 10- 


-4 



'^K = total ph cm ^ s ^ in the line. 
''Equivalent Width 

'^Total phcm-^s"^ in the line from the HETGS observations by Ogle et al. (2000) 



Table 5. Continuum Spectral Fits to the Extended Regions Prom the 3.2 s Frame Time 

Observation 



ModeP 


Region 


Nh (Soft Comp.) 
[xlO^o cin-2] 


kT7 Ti 








Nh (Hard Comp.) 
[xl022cm-2] 


r2 






/ d.o.f. 


I 

II 


NE 

NE 


2.19 (frozen) 
2.19 (frozen) 


9 c+0.2 


(l.4tr2) 


X 10" 
X 10" 


-5 
-5 


2.6 (frozen) 
3.1 (frozen) 


0.32 (frozen) 
0.32 (frozen) 


{6.2tH) X 10- 
(5.5l°i) X 10- 


-6 

-6 


27.3/33 
32.0/33 


I 

II 


SW 

sw 


2.19 (frozen) 
2.19 (frozen) 


q 9+0.1 


(2.ot°o;^) 


X 10" 
X 10" 


-4 
-5 


2.6 (frozen) 
3.1 (frozen) 


0.32 (frozen) 
0.32 (frozen) 


{6.2tli) X 10- 
(5.611 j) X 10- 


-6 
-6 


43.4/48 
54.2/48 



'^Model: I. bremsstrahlung + power law + emission lines; II. power law + power law + emission lines. (See Table 6 for emission line 
parameters.) 
^Unit: keV 

'^Model normalization. For the bremsstrahlung model i^Brom = 3.02 x 10~^^ J nenjdV/ {AnD^), where Ug is the electron density, nj is the 
ion density and D is the distance to the source (all in cgs units). For the power law model, Kpi, = phcm^^ s~^ keV^^ at 1 keV. 



Table 6. Emission Lines Prom the Extended Regions 



ModeP Region Energy Line Observed energy K*' EW 

[keV] [keV] [eV] 



I 


NE 


0.90, 0.91, 0.92 & 1.02 


Ne IX triplet, Ne x Lya 


Qfi+O""' 
^•^"-0.05 




X 10" 


-6 


160. 


I 


NE 


0.74 & 0.77 


VII RRC, VIII Ly/3 


A 79+0.03 
'J- '^-0.04 


(3.9i}:?) 


X 10" 


-6 


101. 


T 
1 




U.OOi, U.OOo & U.O/4 


VII triplet 


U.O»_Q 02 


(q /|+4.2\ 


X iU 


-6 


1 97 


II 


NE 


0.90, 0.91, 0.92 & 1.02 


Nc IX triplet, Ne x Lya 




(4.01?:^) 


X 10" 


-6 


252 


II 


NE 


0.74 & 0.77 


VII RRC, VIII Ly/3 


n 79+0.03 
' ^-0.02 




X 10" 


-6 


172. 


II 


NE 


0.561, 0.568 & 0.574 


VII triplet 


"•OO-0.02 




X 10" 


-5 


194. 


I 


SW 


0.90, 0.91, 0.92 & 1.02 


Ne IX triplet, Ne x Lya 


n QQ+0 02 

u.yo_o.o4 


(4.2tr9) 


X 10" 


-6 


124. 


I 


sw 


0.74 & 0.77 


VII RRC, VIII Ly/J 


0.74100^ 




X 10" 


-6 


86.6 


I 


SW 


0.561, 0.568 & 0.574 


VII triplet 


0.581H? 


(2.61^9) 


X 10" 


-5 


178. 


II 


sw 


0.90, 0.91, 0.92 & 1.02 


Ne IX triplet, Ne x Lya 


q2+° °2 
u.y.^_o.o3 


(5.5ili) 


X 10" 


-6 


221. 


II 


sw 


0.74 & 0.77 


VII RRC, VIII Ly/3 






X 10" 


-6 


180. 


II 


sw 


0.561, 0.568 & 0.574 


VII triplet 


0.58tni 


(o 9+0.5\ 
{'^■'^-0.6) 


X 10" 


-5 


284. 



"^Models: See Table 5. 



= total ph cm ^ s ^ in the line. 
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Fig. 1. — Grey scale representations of Chandra X-ray images of NGC 4151 in the 0.3-2.5 
keV band taken with various frame time. The narrow streak running NE-SW (P.A. 27°) 
across the entire field, and visible in (a) and (b), is not real but a result of the CCD readout 
of the intense nuclear source, (a) 0.1 s frame time. The grey scale is proportional to the 
square root of the X-ray intensity, ranging from cts pixel"^ (white) to 100 cts pixel"^ 
(black), (b) 0.4 s frame time. The grey scale is proportional to the square root of the X-ray 
intensity, ranging from cts pixel"^ (white) to 100 cts pixel"^ (black), (c) 3.2 s frame time, 
the grey scale is proportional to the square root of the X-ray intensity, ranging from cts 
pixel"^ (white) to 71 cts pixel"^ (black). The "hole" at the nucleus in the 3.2 s frame time 
image is not real but caused by pile up. 

Fig. 2. — Contour plot of the 3.2 s frame time Chandra X-ray image of NGC 4151 in the 
0.3-10 keV range. The 11 contour levels range from 4 to 179 cts pixel"^ and are proportional 
to the square root of the X-ray intensity. The hole in the nucleus is not real but a result of 
pile up. 

Fig. 3. — Comparison of model point spread functions (PSFs) with observed radial profiles. 
In each panel, the histogram is the model PSF from the PSF library and is normalized to 
the 0.1 s frame time image peak, which is not piled-up. The points are the observed count 
rates per pixel which were obtained by averaging over an annulus or segments of an annular 
region of radial extent 1". Diamonds represent 0.1 s frame-time data, triangles represent 
0.4 s frame-time data and squares represent 3.2 s frame-time data, (a) The 2-9 keV band. 
The PSF was interpolated to 5.5 keV from the library of PSFs. The observed counts were 
taken from complete, circular annuli and are consistent with the PSF, so the emission in 
this band is unresolved (b) The 0.4 - 2 keV band. The PSF was interpolated to 1.2 keV. 
The observed counts were taken from two pie shaped regions each with opening angle 60° 
along P.A. 144° and P.A. 324°, which are approximately perpendicular to the direction of 
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elongation of the extended emission (see Fig. 1). This diagram shows that the soft emission 
along these directions is marginally resolved or unresolved, (c) Same energy range as (b). 
The counts were taken from a pie shaped region toward the SW (P.A. 233°) with opening 
angle of 120°. (d) Same as (c) but for the NE region (P.A. 54°), with opening angle 120°. 

Fig. 4. — X-ray spectra of the nucleus of NGC 4151, extracted from the 0.1 s and 0.4 s frame- 
time data. The upper panels show the data points with error bars (crosses), with the model 
folded through the instrument response (sohd line). The lower panels show the x residuals 
from this fit. The parameters of these models are listed in Tables 2, 3 and 4. Note that the 
cahbration is uncertain below 0.50 keV and degrades rapidly below 0.45 keV. (a) The 0.1 s 
frame time data with a model comprising two power law continua plus emission lines (Tables 
2 and 3). (b) The 0.1 s frame time data with a model comprising a bremsstrahlung and a 
power-law continuum, plus emission lines (Tables 2 and 3). (c) The 0.4 s frame time data 
with a model comprising a bremsstrahlung and a power-law continuum, plus emission lines 
(Table 4). The channels above 7 keV are ignored in the fit to reduce the effect of "pile- up". 

Fig. 5. — The rectangles defining the SW and NE regions are shown superposed on a grey 
scale representation of the 3.2 s frame-time 0.3 - 10 keV image. Regions within the white 
circle suffer from pile-up. 

Fig. 6. — Hardness ratio (defined as the ratio of the count rates in the 0.8-2 keV to the 
0.4-0.8 keV band) as a function of distance from the nucleus along P.A. 233°. The hardness 
ratio was obtained from the 3.2 s frame time data. These data suffer from pile-up at r < 3". 

Fig. 7. — X-ray spectra of the extended regions of NGC 4151 extracted from the 3.2 s 
frame-time data. The upper panels show the data points with error bars (crosses) and the 
model folded through the instrument response (solid line). The lower panels show the x 
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residuals to this fit. The parameters of the models are hsted in Tables 5 and 6. Note that 
the calibration is uncertain below 0.50 keV and degrades rapidly below 0.45 keV. (a) The 
SW region. The model is a bremsstrahlung plus a power law plus emission lines, (b) The 
SW region. The model is two power laws plus emission lines, (c) The NE region. The model 
is a bremsstrahlung plus a power law plus emission hues, (d) The NE region. The model is 
two power laws plus emission lines. 

Fig. 8. — A superposition of an [OIII] A5007 image (contours, from Perez-Fournon & Wilson 
1990) on the 3.2 s frame time X-ray image. The grey scale is proportional to the square 
root of the X-ray intensity, ranging from 1 cts pixel"^ (white) to 40 cts pixel"^ (black). The 
contour levels for the optical image are drawn at 0.07%, 0.15%, 0.3%, 0.6%, 0.12%, 2.5%, 

5%, 10%, 20% and 80% of the peak. 

Fig. 9. — A superposition of a MERLIN radio continuum image at 21 cm (contours, from 
Mundell et al. 1995) on the 0.1 s frame time grey scale X-ray image. The grey scale is 
proportional to the square root of the X-ray intensity, ranging from cts pixel" ^ (white) to 
a peak intensity of 314 cts pixel" ^ (black). The contour levels of the radio image are drawn 
at 1.5%, 3%, 5%, 12%, 20%, 75% of the peak. 

Fig. 10. — X-ray spectrum of the off nuclear source described in the Appendix. The spectrum 
is extracted from the 3.2 s frame time observation. The upper panel shows the data points 
with error bars (crosses) and the model folded through the instrument response (solid line) . 
The lower panel shows the x residuals to this fit. The model is an absorbed power-law. 
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